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, production of dijets have appeared, we have recalculated these cross sections in next- 

i to-leading order (NLO) of perturbative QCD to see whether they can be interpreted 

s ! ' consistently. The results of these calculations are compared to the data of both collab- 

s ! ' orations. We find that at NLO the cross sections disagree with the data, showing that 

factorization breaking occurs at this order. If direct and resolved contributions are both 
suppressed by the same amount, the global suppression factor depends on the transverse- 
energy cut and is 0.42 for the HI and 0.71 for the ZEUS analysis. However, by suppressing 
only the resolved contribution by a factor of approximately three, also reasonably good 
agreement with all the data is found. The size of the factorization breaking effects for 
resolved photons agrees with absorptive-model predictions. 
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1. Introduction 

It is well known that at high-energy colliders such as the ep collider HERA at 
DESY and the pp collider Tevatron at Fermilab, a large fraction of the observed 
events are diffractive. These events are defined experimentally by the presence of a 
forward-going hadronic system Y with four-momentum py , low mass My (typically 
a proton that remained intact or a proton plus low- lying nucleon resonances) , small 
four-momentum transfer t — [P — Py) 2 , and small longitudinal-momentum transfer 
xp = q(P — py)/{qP) from the incoming proton with four momentum P to the 
central hadronic system X (see Fig. Q] for the case of ep — > eXY). Experimentally, 
a large rapidity gap separates the hadronic system X with invariant mass Mj from 
the final-state system Y with invariant mass My. 

Theoretically, diffractive interactions are described in the framework of Regge 

1 
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Fig. 1. Diffractive scattering process ep — > eXY, where the hadronic systems X and Y are 
separated by the largest rapidity gap in the final state. 



theory W as the exchange of a trajectory with vacuum quantum numbers, the 
pomeron (IP) trajectory. Then the object exchanged between the systems X and 
Y, as indicated in Fig. [T] is the pomeron (or additional lower- lying Regge poles), 
and the upper vertex of the process eJP — ► eX can be interpreted as deep-inelastic 
scattering (DIS) on the pomeron target for the case that the virtuality of the ex- 
changed photon Q 2 = —q 2 is sufficiently large. In analogy to DIS on a proton 
target, ep — > eX, the cross section for the process elP — > eX in the DIS region 
can be expressed as the convolution of partonic cross sections and universal parton 
distribution functions (PDFs) of the pomeron. The partonic cross sections are the 
same as for DIS ep scattering. Usually these pomeron PDFs are multiplied with 
vertex functions for the lower vertex in Fig. [TJ yielding the diffractive parton distri- 
bution functions (DPDFs). The Q 2 -evolution of the DPDFs is calculated with the 
usual DGLAP ^ evolution equations known from ep — > eX DIS. Except for their 
evolution with Q 2 , the DPDFs can not be calculated in the framework of pertur- 
bative QCD and must be determined from experiment. Such DPDFs I 3 | 4 | 5 | 6 | have- 
been obtaine d fr om the HERA inclusive measurements of the diffractive structure 
function F 2 D 12141 defined analogously to the proton structure function Fi . 

Similarly to diffractive DIS, ep — > eXY, where the presence of the large scale Q 
allows for the application of perturbative QCD and X comprises the sum over all 
possible final states, many other processes with a hard scale provided by specific final 
states in the central system X can be predicted using QCD perturbation theory. 
Such processes, usually called hard diffractive processes, are e.g. dijet production 
in diffractive photoproduction (Q 2 ~ 0) and DIS (Q 2 ^ 0), where the large scale 
is given by the jet transverse energy E^ 1 and possibly Q, and diffractive open 
heavy-flavor production, where the large scale is given by the heavy-flavor mass 
and possibly Er and/or Q, in photoproduction or DIS and many more diffractive 
processes induced by pp or pp collisions. The central problem in hard diffraction 
is the problem of QCD factorization, i.e. the question whether diffractive cross 
sections are factorisable into universal DPDFs and partonic cross sections, which are 
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Fig. 2. Diffractive production of dijets with invariant mass M12 in direct (left) and resolved 
(right) photon-pomeron collisions, leading to the production of one or two additional remnant 
jets. 



calculable in perturbative QCD. This question is the subject of the current debate 
in diffractive physics and is of particular interest for the prospects of dis covery of 
new particles such as the Higgs boson in diffractive reactions at the LHC I 7 | 8 | 9 | 

For the inclusive DIS process, factorization has indee d been proven to hold 
and on this basis DPDFs have been extracted at Q 2 ^ EEE from high-precision 
inclusive measurements of the process ep — > eXY using the usual DGLAP evolution 
equations. The proof of the factorization formula, usually referred to as the validity 
of QCD factorization in hard diffraction, also appears to be valid for the production 
of specific final states in DIS, as e.g. the production of jets or heavy-flavor particles, 
and for the direct part of photoproduction {Q 2 ~ 0) or low-Q 2 electroproduction of 
jets^l However, factorization does not hold for hard processes in diffractive hadron- 
hadron scattering. The problem is that soft interactions between the ingoing hadrons 
and/or their remnants occur in both the initial and the final state. This agrees with 
experimental measurements at the TevatronHI Predictions of diffractive dijet cross 
sections for collisions as measured by CDF using DPDFs determined earlier by the 
HI collaboration S2l a t HERA overestimate the measured cross section by up to 
an order of magnitude This large suppression of the CDF cross section can be 
explained by the rcscattcring of the two incoming hadron beams, which, by creating 
additional hadrons, destroy the rapidity gap 

Jet production with real photons involves direct interactions of the photon with 
quarks or gluons originating from the proton or pomeron, respectively, as well as 
resolved photon contributions leading to parton-parton interactions with an addi- 
tional remnant jet coming from the photon as reviewed in El ( see ~pig. [5]). For the 
direct interactions, we expect factorization to be valid as in the case of inclusive 
DIS, as already mentioned, whereas we expect it to fail for the resolved process 
as in hadron-hadron scattering. For this part of photoproduction we would there- 
fore expect a similar suppression factor (sometimes also called rapidity-gap survival 
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probability) due to rescattering effects of the ingoing partons or hadrons. Introduc- 
ing vector- meson dominance photon fluctuations, such a suppressio n b y about a 
factor of three was predicted for resolved photoproduction at HERA 1^-1 

The first measurements of dijet cross sections in diffractive photoproduction 
hav e been presented by the HI collaboration as contributions to two conferences 

The kinematic range for these data were Q 2 < 0.01 GeV 2 , x P < 0.03, E J T etl > 5 
GeV, E^ t2 > 4 GeV and 165 < W < 240 GeV, where jets were identified using the 
inclusive fc-r-cluster algorithm (the definitions of these and the following variables 
will be given in the next section) . The measured cross sections as a function of x° bs 
and Zp s were compared to leading-order (LO) QCD predictions, using the RAP- 
GAP Monte Carlo model 03. For the DPDFs the LO 'HI 2002 fit' was used 
The two cross sections were found to be well described by the predictions in nor- 
malization and shape over the whole range of x° bs and z°p s , showing no breakdown 
of factorization neither in resolved nor in direct photoproduction. In addition, nor- 
malized cross sections as a function of various other variables were compared to the 
predictions with the result that all measured distributions were in good agreement. 

Subsequently we calculated the next-to-leading order (NLO) corrections for the 
cross section of diffractive dijet production usin g the same kinematic cuts and with 
the same DPDFs as in the first HI analysis on the basis of our previous work on 
NLO corrections for inclusive direct^^land resolved^^ dijet photoproduction. While 
at LO good agreement with the HI data was found, consistent with the finding 
in the HI analysis ^^zLit, was found that the NLO corrections increase the cross 
section significantly EQEI] an( j require a suppression factor of the order of R = 0.5. 
Since on theoretical grounds only a suppre ssion of the resolved cross section would 
be acceptable, we demonstrated in | 20 | 21 | ^hat by multiplying the resolved cross 
section with the suppression factor R = 0.34, reasonably good agreement with the 
preliminary HI data could be achieved. This value fo r th e suppression factor 
turned out to be in good agreement with the prediction of ^1 

The first experime ntal data from the ZEUS collaboration were presented at the 
DIS workshop in 2004 ^22l. The dijet cross sections were obtained in the kinematic 
range Q 2 < 1 GeV 2 , x P < 0.035 and E J T etl{2) > 7.5 (6.5) GeV. For these kinematic 
constraints NLO calculations were not available in 2004. So, the measurements were 
compared to LO calculations, unfortunately with previous HI DPDFs ^231 with the 
result, that good agreement in the shape was achieved. But the nor ma lization was 
off by a factor of 0.6, which was attributed to the older DPDF input , so that the 
HI and ZEUS results were consistent with each other. The situation concerning the 
agreement of HI and ZEUS data and the influence of NLO corrections improved 
already considerably in the fall of 2004. These preliminary data from both HERA 
collaborations together with comparisons to NLO calculations based on the DPDF 
fits from were presented at workshops and conferences in the following years. 

In 2006 the HI collaboration published their final DPDF fits from their high- 
precision measurements using the DGLAP evolution equations^. This analysis was 
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based on the larger data sample of the years 1997-2000 as compared to the earlier 
preliminary DPDF sets^. In SI two DPDF sets, the 'HI 2006 fit A' and the 'HI 2006 
fit B' were presented, which both give a good description of the inclusive diffractive 
data. These two sets differ mainly in the gluon density at large fractional parton 
momenta, which is poorly constrained by the inclusive diffractive scattering data, 
since there is no direct coupling of the photon to gluons, so that the gluon density 
is constrained only through the evolution. The gluon density of fit A is peaked at 
the starting scale at large fractional momenta, whereas fit B is flat in this reg ion. 

In 2007 the final publications for diffractive dijet production appeared ^51 The 
comparison between these experimental results and the NLO theory was based on 
the new and final DPDFs from HI 4 . The differential cross sections as measured 
by HI 1251 were compared to NLO predictions obtained with the Frixione program 

interfaced to the 'HI 2006 fit B' DPDFs. The conclusions deduced earlier from 
the comparison with the preliminary data and the preliminary 'HI 2002 fit' ^21 are 
fully confirmed in 1251 with the new DPDFs fits SI. In particular, a global suppression 
is obtained, independent of the DPDFs fits used, i.e. fit A or fit B, by considering 
the ratio of measured dijet cross sections to NLO predictions in photoproduction in 
relation to the same ratio in DIS. In this comparison the value of the suppression 
is 0.5. In addition, by using this overall suppression factor, HI obtained a good 
description of all the measured distributions in the variables Zp s , x° bs , xp, W, 
E J T etl , fj jets , \Ari jets \ and M 12 interfaced with the 'HI 2006 fit B' DPDFs and taking 
into account hadronization corrections!^. Finally, the HI collaboration investigated 
how well the data are describable under the assumption that in the NLO calculation 
the cross section for x° bs > 0.9 is not suppressed. The best agreement in a fit was 
obtained for a suppression factor 0.44 for the NLO calculation with x° bs < 0.9, based 
on fitting the distributions for x° bs 7 W, fp ets and E^?* 1 . In this comparison they 
found disagreement for the largest x° bs -bin and the lowest ff ets (which are related), 
but better agreement in the £^ etl -distribution. In 1251 this leads to the statement, 
that the assumption that the direct cross section obeys factorization is strongly 
disfavored by their analysis. In total, it is obvious that in the final HI analysis I^SI 
a global suppression in diffractive dijet photoproduction is clearly established and 
the model with resolved suppression only is not as well supported by the data. 

Just recently also the ZEUS collaboration presented their final result on diffrac- 
tive dijet photoproduction 123. As in their preliminary analysis, the two jets with 
the highest transverse energies E^ 1 were required to satisfy E^ 11 ^ > 7.5 (6.5) 
GeV, which is higher than in the HI analysis with E 3 T etl{2} > 5 (4) GeV^. ZEUS 
compared their measurements with the NL O predictions for diffractive photopro- 
duction of dijets based on our program 122. Three sets of DPDFs were used, the 
ZEUS LPS fit, determined from a NLO analysis of inclusive diffraction and diffrac- 
tive charm-production data El, and the two HI fits, HI 2006 fit A and fit bSI. The 
NLO results obtained with the two HI fits were scaled down by a factor of 0.87 
4, since the HI measurements used to derive the DPDFs include low-mass proton 
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dissociative processes with My < 1.6 GeV, which increases the photon-diffractive 
cross section by 1.15+q 'Jg as compared to the pure proton final state as corrected 
to in the ZEUS analysis. The comparison of the measured cross sections and the 
theoretical predictions was based on the distributions in the variables y, Mx, xp, 
Zp s , i5jf* , r]f^ and x° bs . The data were reasonably well described in their shape 
as a function of these variables and lay systematically below the predictions. The 
predictions for the three DPDFs differed appreciably. The cross sections for the HI 
2006 fit A (fit B) were the highest (lowest) and the one for the ZEUS LPS fit lay 
between the two, but nearer to the fit A than the fit B predictions. For da jdx obs 
ZEUS also showed the ratio of the data and the NLO predictions using the ZEUS 
LPS fit. It was consistent with a suppression factor of 0.7 independent of x° bs . This 
suppression factor depended on the DPDFs and ranged between 0.6 (HI 2006 fit 
A) and 0.9 (fit B). Taking into account the scale dependence of the theoretical 
predictions the ratio was outside the theoretical uncertainty for the ZEUS LPS fit 
and the HI 2006 fit A, but not for fit B. In their conclusions the authors of the 
ZEUS analysis made the statement that the NLO calculations tend to overesti- 
mate the measured cross section, which would mean that a suppression is present. 
Unfortunately, however, they continued, that, within the large uncertainties of the 
NLO calculations, the data were compatible with the QCD calculations, i.e. with 
no suppression. 

Such a statement clearly contradicts the result of the HI collaboration I^Sl and 
casts doubts on the correctness of the HI analysis. The authors of ^3 attribute this 
discrepancy to the fact that the HI measurements ^Hl we re c arried out in a lower 
E 3 ^ 1 and a higher xp range than those in the ZEUS study ^Zl Besides the different 
E^ 1 and xp regions in and ^3, the two measurements suffer also from different 
experimental cuts of some other variables, which makes it difficult to compare the 
two data sets directly (note also the lower center-of-mass energy for the HI data). In 
addition the comparison with NLO theory in and was done with two different 
programs versus which, however agreed quite well with each other l^S. 

The rather different conclusions concerning factorization breaking in diffract ive 
dijet photoproduction calls for a new comparative study of the two data sets in 1^21 
and ^-U We have therefore performed a new calculation of the NLO cross sections 
on the basis of our earlier work with the new HI 2006 DPDFs and revised 
hadronic corrections as compared to 201 j n ord er t o see whet her we can confirm 
the very different conclusions achieved in the HI ^3 and ZEUS ^3 analyses. In the 
comparison with the new d ata sets we shall follow more or less the same strategy 
as in our earlier workEQEH \\f e first calculate the unsuppressed NLO cross sections 
including an error band based on the scale variation and see how much and in which 
distribution the data points lie inside or outside this error band. Then we determine 
a global suppression factor by fitting the differential cross section da/dE^ 1 at the 
bin with the lowest £5^?* . With this suppression factor we shall compare to the 
differential cross sections of all the other measured variables and look for consistency. 
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In this new comparison between the experimental and the theoretical results we 
shall concentrate on using the HI 2006 fit B SI input, since it leads to smaller cross 
sections than the DPDFs from HI 2006 fit A 4 or the ZEUS LPS DPDF fit 1 

Actually the HI collaboration constructed a third set of DPDFs, which is called 
the 'HI 2007 fit jets'. This fit is obtained through a simultaneous fit to the diffractive 
inclusive and DIS dijet cross sections ^Sl. It is performed under the assumption that 
there is no factorization breaking in the diffractive dijet cross sections. Under this 
assumption, including the diffractive dijet cross sections in the analysis leads to 
additional constraints, mostly on the diffractive gluon distribution. On average the 
'HI 2007 fit jets' is similar to the 'HI 2006 fit B' except for the gluon distribution at 
high momentum fraction and smaller factorization scales. In our analysis we shall 
disregard this new set of DPDFs, since it would be compatible with the factorization 
test of the photoproduction data only if we restricted these tests to the case that 
only the resolved part has this breaking and not the direct part, which has the same 
structure as the DIS dijet cross section. 

In Sec. 2 we shall present the complete list of cuts on the experimental vari- 
ables, give all the input used in the cross section calculations, and present the basic 
formulae, from which the dijet cross sections have been calculated. The comparison 
with the HI ^3 and the ZEUS experimental data is presented and discussed in 
Sec. 3. In this comparison we shall concentrate on the main question, whether there 
is a suppression in the photoproduction data at all. In addition we shall investigate 
also whether a reasonable description of the data is possible with suppression of the 
resolv ed cross section only, as we studied it already in our previous work in 2004 
| 20 | 2 1 | j n g ec ^ we ghgjj fxxiisli with a summary and our conclusions. 

2. Kinematic variables and cross section formulae 
2.1. Kinematic variables and constraints 

The diffractive process ep — * eXY, in which the systems X and Y are separated 
by the largest rapidity gap in the final state, is sketched in Fig. [2j The system X 
contains at least two jets, and the system Y is supposed to be a proton or another 
low-mass baryonic system. Let k and P denote the momenta of the incoming electron 
(or positron) and proton, respectively, and q the momentum of the virtual photon 
7*. Then the usual kinematic variables are 



We denote the four- momenta of the systems X and Y by px and py . The HI data 
^ are described in terms of 



s = 



(fc + P) 2 , Q 2 




(1) 



Ml 



p\ and t = {P-pvY, 

2 , q(P-m 

p Y and xp = — — 




(2) 
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Table 1. Kinematic cuts applied in the HI 
analysis of diffractive dijet photo-production. 



165 GeV < 


W 


< 


242 GeV 




Q 2 


< 


0.01 GeV 2 




11/ rp 


> 


5 GeV 






> 


4 GeV 


-1 < 


jetl,2 
^lab 


< 


2 




Xp 


< 


0.03 




M Y 


< 


1.6 GeV 




-t 


< 


1 GeV 2 



where Mx and My are the invariant masses of the systems X and Y, t is the 
squared four-momentum transfer of the incoming proton and the system Y, and 
xp is the momentum fraction of the proton beam transferred to the system X . 

The exchange between the systems X and Y is supposed to be the pomeron 
IP or any other Regge pole, which couples to the proton and the system Y with 
four-momentum P — py . In this work we will neglect Reggeon exchanges, which 
contribute only at large xp . The pomeron is resolved into partons (quarks or gluons) 
with four- momentum v. In the same way the virtual photon can resolve into partons 
with four-momentum u, which is equal to q for the direct process. With these two 
momenta u and v we define 

Pu , qv 
Pq q(P-py) 

x~f is the longitudinal-momentum fraction carried by the partons coming from the 
photon, and zp is the corresponding quantity carried by the partons of the pomeron 
etc., i.e. the diffractive exchange. For the direct process we have = 1. The 
final state, produced by the ingoing momenta u and v, has the invariant mass 
M\2 = y (u + v) 2 , which is equal to the invariant dijet mass in the case that no more 
than two hard jets are produced, q — u and P — py — v are the four-momenta of the 
remnant jets produced at the photon and pomeron side. The regions of the kinem atic 
variables, in which the cross section has been measured by the HI collaboration 
are given in Tab. [lj whereas the corresponding regions for the ZEUS analysis 
are given in Tab. 2. In each case, we have evaluated the theoretical cross sections 
with the corresponding constraints. 

The upper limit of x p is kept small in order for the pomeron exchange to be 
dominant. In the experimental analysis as well as in the NLO calculations, jets are 
defined with the inclusive fc^-cluster algorithm with a distance parameter d = 1 
1221 in the laboratory frame. At least two jets are required with transverse energies 
E 3 T etl > 5 (7.5) GeV and E J T et2 > 4 (6.5) GeV. They are the leading and subleading 
jets with -1 < rf^l 1 ' 2 < 2 (-1.5 < rf^' 2 < 1.5) for HI (ZEUS). The lower limits of 
the jet -Et's are asymmetric in order to avoid infrared sen sitiv ity in the computation 
of the NLO cross sections, which are integrated over Et ^1 

In the experimental analysis the variable y is deduced from the energy E' e of the 
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Table 2. Kinematic cuts applied in the ZEUS 
analysis of diffractive dijet photoproduction. 



0.2 


< 


y 


< 


0.85 






Q 2 


< 


1 GeV 2 






E{ etl 


> 


7.5 GeV 








> 


6.5 GeV 


-1.5 


< 


jetl,2 


< 


1.5 






Xp 


< 


0.025 






-t 


< 


5 GeV 2 



scattered electron, y = 1 — E' e /E e . Furthermore, sy — W 2 = (q + P) 2 = (px + py) 2 . 
The range of W given in Tab. 1 corresponds to the y range 0.3 < y < 0.65. xp is 
reconstructed according to 

where E p is the proton beam energy and the sum runs over all particles (jets) in the 
X-system. The variables M12, x-y, and zp are determined only from the kinematic 
variables of the two hard leading jets with four- momenta p> etl and p> et2 . So, 

^i2 = (P ictl +P ict2 ) 2 , (5) 
where additional jets are not taken into account. In the same way we have 



= ^^-p,) and = ^ E+P *>. (6 ) 
2yE e F 2xpE p 



The sum over jets runs only over the variables of the two leading jets. These defini- 
tions for x 7 and zp are not the same as the definitions given earlier, where also the 
remnant jets and any additional hard jets are taken into account in the final state. 
In the same way M x can be estimated by M\ — M 2 2 / (zp^ s x° hB ). The dijet system 
is characterized by the transverse energies E 3 ^ tl and E^ t2 and the rapidities in the 
laboratory system 77^ and r^ 2 . The differential cross sections are measured and 
calculated as functions of the transverse energy Eip of the leading jet, the average 
rapidity fj jets = {rf^f + rf { *f)/2, and the jet separation IAt^I = \rf^f - rf { *f\, 
which is related to the scattering angle in the center-of-mass system of the two jets. 



2.2. Diffractive parton distributions 

The diffractive PDFs are obtained from an analysis of the diffractive process 
ep — > eXY, which is illustrated in Fig. [TJ where now Q 2 is large and the state 
X consists of all possible final states, which are summed. The cross section for this 
diffractive DIS process depends in general on five independent variables (azimuthal 
angle dependence neglected): Q 2 , x (or xp, My, and t. These variables are 
defined as before, and x = Q 2 /(2Pq) = Q 2 /(Q 2 + W 2 ) = x P (3. The system Y is 
not measured, and the results are integrated over — t < 1 GeV 2 and My < 1-6 GeV 
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as in the photoproduction case. The measured cross section is expressed in terms 
3 c 

d 3 cr D ina 2 



of a reduced diffractive cross section oy 5 ^ defined through 



dxpdxdQ 2 xQ 4 



*0 a^\x P ,x,Q 



l-y+^)a^(x P ,x,Q 2 ) (7) 



and is related to the diffractive structure functions F 2 D ' 3 ' and F^ 3 ^ by 

D(3) _ pD(3) _ V 2 F D(3) fgl 

°r ~*2 l + (l-y)2^ " ^ 

y is defined as before, and is the longitudinal diffractive structure function. 

The proof of Collins that QCD factorization is applicable to diffractive DIS, 
has the consequence that the DIS cross section for j*p — » XY can be written as a 
convolution of a partonic cross section , which is calculable as an expansion in 
the strong coupling constant a s , with diffractive PDFs yielding the probability 
distribution for a parton a in the proton under the constraint that the proton 
undergoes a scattering with a particular value for the squared momentum transfer 
t and xp. Then the cross section for j*p — > XY is 

^ = E/7w(^Q 2 > ^feQ 2 ;^,t). (9) 

This formula is valid for sufficiently large Q 2 and fixed xp and t. The parton cross 
sections are the same as those for inclusive DIS. The diffractive PDFs are non- 
perturbative objects. Only their Q 2 -evolution can be predicted with the well-known 

l9l 

DGLAP evolution equations , which we shall use in NLO. 

Usually for f®(x,Q 2 ;xp,t) an additional assumption is made, namely that it 
can be written as a product of two factors, fp/ p (xp,t) and f a /p(P, Q 2 ), 

fa{ x iQ 2 'i x P^) = fp/p(xp,t)f a /p{f3 = x/xp,Q 2 ). (10) 

fp/p( x Pi t) is the pomeron flux factor. It gives the probability that a pomeron with 
variables xp and t couples to the proton. Its shape is controlled by Regge asymp- 
totics and is in principle measurable by soft processes under the condition that 
they can be fully described by single-pomeron exchange. This Regge factorization 
formula represents the resolved pomeron model, in which the diffractive exchange, 
i.e. the pomeron, can be considered as a quasi-real particle with a partonic struc- 
ture given by PDFs f a /p{P,Q 2 )- P is the longitudinal momentum fraction of the 
pomeron carried by the emitted parton a in the pomeron. The important point is 
that the dependence of J® on the four variables x,Q 2 ,xp and t factorizes into two 
functions fp/ p and f a /p, which each depend only on two variables. 

Since the value of t could not be fixed in the diffractive DIS measu remen ts, it is 
integrated over with t in the region t cut <t< i min . Therefore we have^ES! 

f(xp)= / dtfp /p (xp,t), (11) 

"'tout 
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w here t cut = — 1 GeV 2 and i m i n is the minimum kinematically allowed value of \t\. 
In ™ the pomcron flux factor is assumed to have the form 

f P/p (x P) t) = x l p 2aip{t) exp(B P t). (12) 

ap(t) is the pomeron trajectory, oip(t) = ap(0) + a'pt, assumed to be linear in t. 
The values of Bp. ap(0) and a'p are taken from SI and have the values Bp = 5.5 
GeV" 2 , a P {0) = 1.118 (fit A), a P {0) = 1.111 (fit B) and a' P = 0.06 GeV" 2 . 
Usually fp/ p {xp,t) as written in Eq. (jT2)) has in addition to the dependence on xp 
and t a normalization factor N, which can be inferred from the asymptotic behavior 
of a tot for pp and pp scattering. Since it is unclear whether these soft diffractive cross 
sections are dominated by a single pomeron exchange, it is better to include N into 
the pomeron PDFs f a /p and fix it from the diffractive DIS data SI. The diffractive 
DIS cross section cy* is measured in the kinematic range 3.5 < Q 2 < 1600 GeV 2 , 
0.01 < [3 < 0.9 and 10~ 4 < x P < 0.05. 

The pomeron couples to quarks in terms of a light flavor singlet S(zjp) = u(zp) + 
d{zp) + s{zp) + u{zp) + d(zp) + s(zp) and to gluons in terms of g{zp), which are 
parameterized at the starting scale Qg = 1.75 GeV 2 (fit A) and 2.5 GeV 2 (fit B). zp 
is the momentum fraction entering the hard subprocess, so that for the LO process 
zp = {3, and in NLO < Zp < 1. These PDFs of the pomeron are parameterized 
by a particular form in terms of the usual power ansatz as given in ^. Charm quarks 
and bottom quarks couple differently from the light quarks by including the finite 
charm mass m c = 1.4 GeV and bottom mass rrn, = 4.5 GeV in the massive quark 
scheme and describing the coupling to photons via the photon-gluon fusion type 
process up to order a 2 . For the pomeron PDFs, we used a two-dimensional fit in 
the variables zp and Q 2 and then inserted the interpolated result in the cross section 
formula. 



2.3. Cross section formula 

Under the assumption that the cross section can be calculated from the well-known 
formulae for jet production in low-Q 2 ep collisions, the cross section for the reaction 
e + p — > e + 2 jets + X' + Y is computed from the following basic formula: 

da D (ep + 2 jets + X' + Y) = 

y I dt dxp / dzp / dy dx 1 

a b Jt^t Jx>™ Jo /y min JQ 

f 1 /e{y)fa/ 1 (x 1 ,M 2 )fp/ p (xp,t)f b/ p(zp, Mp) 

da^(ab ->jets). (13) 

y, Xj and zp denote the longitudinal momentum fractions of the photon in the 
electron, the parton a in the photon, and the parton b in the pomeron. M 1 and Mp 
are the factorization scales at the respective vertices, and da^ n \ab — > jets) is the 
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cross section for the production of an n-parton final state from two initial partons 
a and b. It is calculated in NLO, as are the PDFs of the photon and the pomeron. 

The function / 7 / e (y), which describes the virtual photon spectrum, is assumed 
to be given by the well-known Weizsacker- Williams approximation, 



f,/e(y) = J 



1 + (1 - V? ln Qmax(l ~ V) 



y m^y 2 



2m\y 



l-y 1 



m e2/ 2 Qmax 



(14) 



Usually, only the dominant leading logarithmic contribution is considered. We have 
added the second non-logarithmic term as evaluated in^H. Q^ax — 0-01 (1) GeV 2 
for the HI (ZEUS) cross sections calculated in this work. 

The formula for the cross section da D can be used for the resolved as well as for 
the direct process. For the latter, the parton a is the photon and f 1 / 1 (x 1 , M 2 ) = 
<5(1 — Xj), which does not depend on M 7 . As is well known, the distinction between 
direct and resolved photon processes is meaningful only in LO of perturbation the- 
ory. In NLO, collinear singularities arise from the photon initial state, that must be 
absorbed into the photon PDFs and produce a factorization scheme dependence as 
in the proton and pomeron cases. The separation between the direct and resolved 
processes is an artifact of finite order perturbation theory and depends in NLO on 
the factorization scheme and scale M 7 . The sum of both parts is the only physically 
relevant quantity, which is approximately independent of the factorization scale M 7 
due to the compensation of the scale dependence between the NLO direct and the 
LO resolved contribution. For the resolved process, PDFs of the photon are needed, 
for which we choose the NLO versions of GRV transformed to the MS scheme. 



3. Results 

3.1. Comparison with HI data 

In this section, we present the comparison of the various theoretical predictions in 
NLO with the experimental data from the HI collaboration^^. The corresponding 
kinematic cuts are given in Tab. 1. Before we confront the calculated cross sec- 
tions with the experimental data, we correct them for hadronization effects. The 
hadronization corrections are calculated by means of the LO RAPGAP Monte Carlo 
generator^^. The factors for the transformation from jets made up of stable hadrons 
to parton jets were supplied by the HI collaboration ° . Most of our calculations are 
done with the 'HI 2006 fit B DPDFs since they give they smaller diffractive dijet 
cross sections as compared to the 'HI 2006 fit A'. These DPDF fits are based on 
Uj = 3 massless flavors. The production of charm and bottom quarks was treated 
there in the Fixed-Flavor Number Scheme (FFNS) in NLO with non-zero charm 
and bottom quark mass. Instead of this extra treatment of the charm and bottom 
contribution in the pomeron we added a charm PDF in the pomeron as obtained 
in the 'HI 2002 fit' , where the charm quark was also considered to be massless. 
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The bottom contribution was neglected. This assumption simplifies the calculations 
considerably. Since the charm contribution from the pomeron is small, this should 
be a good approximation. We then take nf = 4 with A^ = 0.347 GeV, which 

corresponds to the value used in the DPDFs 'HI 2006 fit A' and 'HI 2006 fit B'Sl. 

As it is clear from the discussion of the various preliminary analyses of the HI and 
ZEUS collaborations, there are two questions which we would like to answer from 
the comparison with the recent HI and the ZEUS data. The first question is whether 
a suppression, which differs substantially from one, is needed to describe the data. 
The second question is whether the data are also consistent with a suppression factor 
applied to the resolved cross section only. To give an answer to these two questions 
we calculated first the cross sections with no suppression factor (R — 1 in the 
following figures) with a theoretical error obtained from varying the common scale 
of renormalization and factorization by factors of 0.5 and 2 around the central scale 
(highest E^). In a second step we show the results for the same differential cross 
sections with a global suppression factor, adjusted to da/dE^ 1 in the smallest 
_E^ etl -bin. As in the experimental analysis ^Sl, we consider the differential cross 
sections in the variables if, zp , log 1Q (x F ), E 3 T etl , M 12 , f] jets , \A?f ets \ and W. 

The unsuppressed (R = 1) cross sections da/dx° bs , da/dzp s , da/d\og w (x]p), 
d<r/dE J T et \ da/dM 12 , da/di-f ets , dcr/d\A^ ets \ and da/dW {ff ets = (tjf^) in ESI) 
with their scale variation are shown in Fig. 3a-h. In these figures we also plotted 
the experimental data with their errors. Except for two points (largest Zp s and 
largest £^ etl -bin) all other experimental points lie, including their errors, outside 
the theoretical error band. This comparison clearly tell us, that an unsuppressed 
cross section is in disagreement with the data. It is clear, that with the DPDFs 'HI 
2006 fit A' cross section this conclusion would be even stronger, since with these 
DPDFs the unsuppressed cross sections are even larger. That da/dzp s overlaps in 
the largest bin with the lower limit of the prediction for R = 1 (see Fig. 3b) can be 
explained with the fact that the gluon DPDF in the 'HI 2006 fit B' is not very well 
constrained for large [3 and might be larger there. 

If we now determine the suppression factor from fitting the lowest £^ etl -bin 
experimental cross section we obtain R = 0.42 ± 0.06. The indicated error corre- 
sponds to the experimental uncertainty, while we show in the figures explicitly the 
theoretical uncertainty. With this suppression factor we have calculated the eight 
distributions including their theoretical errors and compare with the experimental 
data including their errors. The results of this comparison is shown also in Figs. 
3a-h. With the exception of Figs. 3d and 3h, where the comparisons of da/dE^ n 
and da/dMi2 are shown, all other plots are such that the data points lie inside the 
error band based on the scale variation. Most of the data points even agree with 
the R = 0.42 predictions inside the much smaller experimental errors. In da/dE 3 T etl 
(see Fig. 3d) the predictions for the second and third bin lie outside the data points 
with their errors. For R = 1 and R = 0.42 these cross sections falls off stronger 
with increasing E°^ tl than the data, the normalization being of course about two 
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Fig. 3. Differential cross sections for diffractive dijet photoproduction as measured by HI and 
compared to NLO QCD without (R = 1) and with (R = 0.42) global suppression (color online). 



times larger for R = 1. In particular, the third data point agrees with the R = 1 
prediction. This means that the suppression decreases with increasing E^ a . Such 
a behavior points in the direction that a suppression of the resolved cross section 
only would give better agreement with the data, as we shall see below. The same 
observations can be made by looking at dajdMyi in Fig. 3e. The survival proba- 
bility R = 0.42 ± 0.06 agrees with the result in 1221, which quotes R = 0.50 ± 0.10, 
determined by a fit to the double ratio of measured to predicted cross section in 
photoproduction by the corresponding ratio in DIS given as a function of W. In this 
double ratio many experimental errors and theoretical scale errors cancel to a large 
extent. This double ratio is also insensitive to the detailed shape of the diffractive 
gluon density. From our comparison we conclude that the HI data sho w a gl obal 
suppression of the order of two in complete agreement with the results | 20 | 21 | ari d 
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^3 based on earlier preliminary an( j final HI data 

Next we want to answer the second question, whether the data could be con- 
sistent with a suppression of the resolved component only. For this purpose we 
have calculated the cross sections in two versions: (i) suppression of the resolved 
cross section and (ii) suppression of the resolved cross section plus that part of 
the NLO direct part which depends on the factorization scale at the photon ver- 
tex and ther eby eliminates the My-dependence in the combined direct and resolved 
cross section Of course, the needed suppression factors for the two versions will 
be different. We determine the suppression factors by fitting again the measured 
da/dE^ n for the lowest i^ etl -bin (see Fig. 4d). Then, the suppression factor for 
version (i) is R = 0.31 (denoted res in the figures), and for version (ii) it is R = 0.29 
(denoted res+dir-IS). The comparison with the HI data of da/dx obs 7 da/dz$ s , 
da/dlog w (x P ), da/dE 3 T et \da/dM 12 , da/drj 3ets , da/d\A^ ets \ and da/dW is shown 
in Figs. 4a-h, where we also have plotted the prediction for the global suppression 
(direct and resolved) with R = 0.42, already shown in Figs. 3a-h. Looking at Figs. 
4a-h we can distinguish three groups of results from the comparison with the data. 
In the first group, the cross sections as functions of Zp s , log 10 (a;_p), Mi 2 , \Arj 3ets \ 
and W, the agreement with the global suppression (R = 0.42) and the resolved 
suppression (R = 0.31 and R = 0.29) is comparable. In the second group, which 
consists just of da / 1 dEi^ tl , the agreement is better for the resolved suppression only. 
In the third group, dajdx obs and da/dr) jets , the agreement with the resolved sup- 
pression is worse than with the global suppression. In particular, for da/dx° bs , which 
is usually considered as the characteristic distribution for distinguishing global ver- 
sus resolved suppression, the agreement with resolved suppression does not improve. 
Unfortunately , th is cross section has the largest hadronic corrections of the order 
of (25 — 30)% 1251. Here, the bins with largest x° bs are particularly sensitive to the 
hadronic corrections and possible migrations of the data between the two bins. If 
we average the cross sections for these two bins, the agreement with the data point 
becomes much better. We also notice, that the predictions for the two suppression 
modes (i) and (ii) are almost the same. The only exception are the cross sections for 
the largest ir° fcs -bin (see Fig. 4a). In Figs. 4a-h the theoretical errors coming from 
the scale uncertainty are not shown. If they are taken into account, the difference 
between experimental data and theory in Figs. 4a and 4f is much less dramatic. 
On the other hand, for the cross section da /dE^ tl the agreement improves con- 
siderably with the suppression of the resolved part only (note the logarithmic scale 
in Fig. 4d). Here, of course, we must admit that the suppression factor could be 
i?T-dependent, although we do not know of any mechanism, which could cause such 
a E^-dependence of the suppression. We remark that this .©^'-dependence of the 
global suppression is also visible in the HI analysis of ^3. 

We also checked for two distributions whether the predictions for resolved sup- 
pression depend on the chosen diffractive PDFs. For this purpose we have calculated 
for the two cases da/dzp 8 and da/dE J ^ tl the cross sections with the 'HI 2006 fit 
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Fig. 4. Differential cross sections for diffractive dijet photoproduction as measured by HI and 
compared to NLO QCD with global, resolved, and resolved/direct-IS suppression. 



A' partem distributions^. The results are compared in Figs. 5a and b to the results 
with the 'HI 2006 fit B' and the experimental data. Of course, since the 'HI 2006 
fit A' PDFs have a larger gluon component at large /?, the cross sections are larger 
and therefore need a larger suppression factor R = 0.22. From Figs. 5a and b we 
conclude that there is no appreciable dependence on the chosen DPDFs. Note that 
in Fig. 5b the cross section for the smallest .E^f^-bin has been fitted to determine 
the suppression factor. In total, we are tempted to conclude from the comparisons 
in Figs. 4a-h that the predictions with a resolved-only (or resolved+direct-IS) sup- 
pression are consistent with the HI data ^4 The only exceptions are two bins in 
the x° bs and one bin in the ^^-distribution. 
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Fig. 5. Differential cross sections for diffractive dijet photoproduction as measured by HI and 
compared to NLO QCD with resolved suppression and two different DPDFs. 



3.2. Comparison with ZEUS data 

In this subsection we shall compare our predi ctio ns with the final analysis of the 
ZEUS data, which was published just recently 123. The kinematic cuts are given in 
Tab. 2. There are the following differences to the HI cuts in Tab. 1: First the upper 
cut on Q 2 is larger. Second there is a larger range in the variable y and the upper 
cut on xp is slightly smaller. The most important change is the larger IS^f* 1 ' 2 ' cut, 
namely E 3 ^ 1 ^ > 7.5 (6.5) GeV, which leads to smaller cross sections. Also the cut 
on \t\ is different. The different cuts on Q 2 and |i| have little influence. For example, 
the larger |t|-cut in Tab. 2 as compared to Tab. 1 increases the cross section only 
by 0.2%. The constraint on My is not explicitly given in the ZEUS publication 123 
They give the cross section for the case that the diffractive final Y state consists only 
of the proton. For this they correct their measured cross section by subtracting in all 
bins the estimated contribution of a proton-dissociative background of 16%. When 
comparing to the theoretical predictions with the DPDFs from the HI 2006 fits, 
they multiply the theoretical cross section with a (slightly different) factor of 0.87 
in order to correct for the proton-dissociative contributions, which are contained 
in these DPDFs by requiring My < 1.6 GeV. We do not follow this procedure. 
Instead we leave the theoretical cross sections unchanged, i.e. they contain a proton- 
dissociative contribution with My < 1.6 GeV and multiply the ZEUS cross sections 
by 1.15 to include the proton-dissociative contribution. In this comparison we shall 
follow the same strategy as before. Before we do this, we correct our theoretical 
prediction by the hadronization corrections reported in 123. We first compare to the 
predictions with no suppression (i? = 1) and then determine a suppression factor by 
fitting da/dE^ n at the smallest iXf' -bin. Then we compare to the cross sections 
as a function of the seven observables x° bs , Zp s , xp, £|f* , y, Mr and r) jetl instead 
of the eight variables in the HI analysis. The distribution in y is equivalent to the 
TU-distribution in 1213. The ZEUS collaboration have also published experimental 
measurements in the two regions x~ t < (>) 0.75, which do however not consider 
here due to space limitations. 
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Fig. 6. Differential cross sections for diffractive dijet photoproduction as measured by ZEUS and 
compared to NLO QCD without (Ft = 1) and with (Ft = 0.71) global suppression (color online). 



The theoretical predictions for these differential cross sections with no suppres- 
sion factor (R — 1) are shown in Figs. 6a-g, together with their scale errors and 
compared to the ZEUS data points. Except for the x° bs and £yf distributions, 
most of the data points lie outside the theoretical error bands for R = 1. In par- 
ticular, in Figs. 6b, c, e, f and g, 2, 3, 4, 4 and 5 points lie outside. This means 
that most of the data points disagree with the unsuppressed prediction. Next, we 
determine the suppression factor from the measured da/dE^ 1 at the lowest E^ 11 - 
bin, 7.5 GeV < E 3 T etl < 9.5 GeV, and obtain R = 0.71 ± 0.06. The indicated error 
corresponds again to the experimental uncertainty, while we show in the figures 
explicitly the theoretical uncertainty. This means that the suppression factor from 
the ZEUS data is larger than the one obtained from the analysis of the HI data, 
which is actually consistent with the da/dEi^ 11 for the second bin in Fig. 3b. Here 
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the cross section is approximately larger by a factor of 1.8 than the prediction with 
R = 0.41. If we now check how the predictions for R = 0.71 compare to the data 
points inside the theoretical errors, we observe from Figs. 6a-g that with the ex- 
ception of da/dzp s (largest bin), the data points agree with the predictions inside 
the theoretical error band. This is quite consistent with the HI analysis, discussed 
in the previous subsection, and leads to the conclusion that also the ZEUS data 
agree much better with the suppressed predictions than with the unsuppressed one. 
In particular, the global suppression factor approximately agrees with the global 
suppression factor, which one would expect from the analysis of the HI data at the 
second smallest E^ n -bm. 

Similarly as in the previous section we compared the ZEUS data also with the 
assumption that the suppression results only from the resolved cross section. Here we 
consider again the two versions: (i) only resolved suppression (res) and (ii) resolved 
plus direct suppression of the initial-state singular part (res+dir-IS). For these two 
models we obtain the suppression factors R — 0.53 and R = 0.45, respectively, 
where these suppression factors are again obtained by fitting the data point at the 
first bin of da/dE^* . The comparison to the global suppression with R = 0.71 and 
to the data is shown in Figs. 7a-g. In general, we observe that the difference between 
global suppression and resolved suppression is not large, i.e. the data points agree 
with the resolved suppression almost as well as with the global suppression. 

In Figs. 8a and b the difference between 'HI 2006 fit B' and 'HI 2006 fit A' 
is shown again for the case of the resolved suppression. In both figures we observe 
that the fit A suppression with the suppression factor R = 0.27 agrees better with 
the data than with the factor R = 0.53 for the fit B suppression. In particular, for 
da/dE 3 ^ 1 the agreement with the three data points is perfect (note the logarithmic 
scale). 

4. Conclusion 

In summary, we have revisited the final HI and ZEUS data on the diffractive photo- 
production of dijets at HERA. We focused on the question if the two data sets, taken 
with different ep center-of-mass energies and kinematic cuts (in particular on the 
jet transverse energies), could be consistently interpreted within QCD factorization, 
employing universal parton densities in the diffractive exchange and process-specific 
hard partonic cross sections evaluated at NLO, or showed rather evidence of factor- 
ization breaking in the direct and/or resolved photon channels. 

First, we found that even with the most optimistic (and likely realistic) parton 
density set 'HI 2006 fit B', both the HI and ZEUS data sets were overestimated 
by the unsuppressed NLO predictions and better described by global suppression 
factors of 0.42 ± 0.06 and 0.71 ± 0.06, respectively. These factors were obtained by 
fitting our NLO predictions in both cases to the lowest (and dominant) £^ e *-bin 
and are in agreement with the global suppression factor of 0.50 ± 0.10 found by the 
HI collaboration in a fit to all of their data points, but at variance with the final 
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Fig. 7. Differential cross sections for diffractive dijet photoproduction as measured by ZEUS and 
compared to NLO QCD with global, resolved, and resolved/direct-IS suppression. 



conclusion in the publication by the ZEUS collaboration. 

Second, we demonstrated that the HI (ZEUS) data (in particular the Et- 
distributions and somewhat less the rr 7 -distributions, which are unfortunately sub- 
ject to large hadronization uncertainties) can be described almost equally well by 
applying a suppression factor of about one-third (one-half) to the resolved-only 
contribution. We showed that this could be consistently done by suppressing also 
the direct intial-state singular part without a big impact on the suppression factor 
with the added advantage of preserved factorization-scale invariance. Alternatively, 
we admitted for the possibility that a global suppression factor might be in- 
dependent, although a theoretical motivation is only known for the first scenario 
and the suppression factor obtained is in good agreement with absorptive-model 
predictions. 
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Fig. 8. Differential cross sections for diffractive dijet photoproduction as measured by ZEUS and 
compared to NLO QCD with resolved suppression and two different DPDFs. 



Finally, we showed that our conclusions, while the numerical values of the sup- 
pression factors may change to some extent, are not qualitatively altered, when a 
different set of diffractive parton densities (e.g. 'HI 2006 fit A') is employed. The 
same observation should hold for the very recent 'HI 2007 fit jets', which is very 
similar to the 'HI 2006 fit B'. 

While the epoch of HERA experiments has now ended and an International 
Linear Collider may not be built in the near future, it will be very interesting to 
investigate diffractive physics at the LHC. As stated above, the search for Higgs 
bosons may benefit in an important way from the diffractive production channels, 
and this depends crucially on an excellent understanding of the QCD backgrounds. 
Proton-proton and heavy-ion collisions at the LHC may even be a source of high- 
energy photon collisions, and this may open up a whole new field of investigation 
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for diffractive dijet photoproduction 1 ^. 
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